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Abstract 

Methods of classical perturbation theory developed for small 
perturbations are extended to slowly (or adiabatically) perturbed 
systems, with slow dependence either on time or on dynamical vari- 
ables. Specifically, the extension is performed for the canonical 
perturbation theory of Poincare and Von Zeipel, for the Krylov- 
Bogoliubov-Kruskal method of eliminating angle variables, for the 
general form of direct near-identity canonical transformations and 
for two of its realizations, based on the "conventional" generating 
function and on the Lie transform. In addition, the concepts of slow 
(or adiabatic) perturbations and of an implicit "small parameter" £ 
are clarified, as is the distinction between two alternative defi- 
nitions of adiabatic invariance, and as an example the solution of 
the slowly perturbed harmonic oscillator up to and including 0( £/*) 
is derived. 
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INTRODUCTION 




Perhaps the most widely studied perturbation problem in classical mechanics 
is that of perturbed periodic motion . If a motion is given that is soluble and 
periodic, the problem may be concisely defined as the derivation of an approxi- 
mate solution for a motion that is slightly different. 


This "slight change" applied to the motion is termed the perturbation and 

(or "adiabatic") 

it usually belongs to one of two types: " small " perturbations and " slow" ^ ones. 
The difference between the two is best explained by assuming that the motion can 
be described by a Hamiltonian, although this condition is not essential. In a 
slightly perturbed motion the Hamiltonian may then be written 


H = H^ + £H^ + £ 2 + ... 


( 1 ) 


where £ « 1 is a small numerical parameter characterizing the magnitude 
of the perturbation and where the limit L 0 corresponds to the unper- 
turbed motion. A typical example v;;>uld be the motion of a planet around the 

(o) 

sun as perturbed by the planet Jupiter. In that case H describes 'die 

planet's Keplerian motion in the sun's gravity field while H^^ describes the 

lowest order of the perturbation induced by Jupiter. The zero-order Hamiltonian 

is then proportional to the solar mass mg while is proportional 

to the mass m. of Jupiter: the ratio of the two terms will be of the order 

(m./m ) (about 10 - ^) and this dimensionless quantity provides a natural 
«} 8 

choice for £ . 


To illustrate a slow perturbation, consider a Hamiltonian that is slowly 
dependent on the time t (slow dependence may also involve canonical, variables): 


i «_ 


H = H(p, q, t) 


( 2 ) 



K 


Then the dependence is said to be slow If the terms produced by the opera- 
tion "d/ut are by an order in £ smaller than the terms from which they 
are derived, e.g. 

n>H/~6t a 0(£H) (3) 

The preceding equation is not quite precise, since it implies that £. has 
the dimension of t~^ . In fact, one always requires some natural time period 
T against which the rapidity of the time variation may be gauged, this usually 
being the period of the unperturbed system. With this taken into account, (3) 
becomes 

n>H/t)t = 0 (€-H/t) ( 4) 

and 6 is clearly dimensionless. 

In either type of problem there generally exists a steadily increasing 
"angle variable" appear! g in the argument of sines and cosines, describing 
the nearly- periodic part of the motion. 0: a way of "solving" the problem then 
involves finding a transformation to new variables, such that the angle 
variable is eliminated from the equations of motion. If the system also 
possesses a Hamiltonian H , the absence of the angle variable from H implies 
that its conjugate "action variable" is a constant of the motion, and this 
eliminates an additional variable from consideration. In slowly perturbed 
systems, such constants are called adiaba tic invarian t s ^ ^ . In slightly perturbed 

systems, no generally accepted name exists (G. Contopoulos, who investigated 

( 2 ) 

the relation between the two types of constants has termed them "third 
integrals") but they are well-known in celestial mechanics and may be derived 


in a variety of ways 


Tho purpose of this work is to show how the standard methods of celestial 
mechanics, designed to handle small perturbations, may be modified to deal 
with slow perturbation and lead to the derivation of adiabatic invariants. 

Two methods will be considered here: the PoincarA-Von Zeipel method''’'"' ' * 

(7) 

for solving the Kami 1 t on-J acobi equation and the Krylov-Bogoliubov procedure 
(or the related method of Kruakal) . In addition, it will be shown that the 
direct form of near-identity canonical transformations can also be adapted 
to cases in which some variables are slow. 

EXPLICIT AND IMPLICIT 6 

In the example of perturbed planetary motion the small parameter £ can be 
given an explicit numerical value. In problems of slowly perturbed motion this 
is often difficult to do and one may then speak of an implicit £ * 

As the arch type of a slowly perturbed system, consider the "pulled-up 
pendulum. »^5)(l4) J a gj^p^e pendulum is suspended from a hole in the ceiling 

and its suspension string is pulled up (or released) at a slow, though not 

to 

necessarily constant rate. Obviously, the angular frequency^of the pendulum 
will vary and, since work is being done against the centrifugal force of the 
oscillation, so will its energy E . However, as long as the rate at which 
the string is withdrawn is sufficiently slow (and does not resonate with the 
oscillation of the pendulum) an adiabatic invariant may be found, equaling 
E/<o in the lowest order. 

Two points should be noted here. First, the perturbation need not be 
small: by the time the withdrawal is complete, the length of the pendulum 
may well have changed by a large factor. Secondly, while one can devise an 
explicit 6 for the problem — e.g. £ ** , where t is the time in 


which the length of the pendulum la reduced to l/e of its value, at the 

given (instantaneous) rate — its value nowhere enters the calculation. 

A more complicated example is provided by the motion of a charged 

particle in a slightly inhomogeneous magnetic field B • Here "slightly" 

means that the derivatives "D B./"6 x . are all of order £ smaller than 

i j 

the components of the field intensity and its magnitude B • Thus the slow- 
ness is in the dependence on spatial coordinates and a scale-length for 
gauging it is provided by the gyration radius $ , giving, in analogy to 

eq. ( 4 ) 

^ x. * 0(E b/? ) (5) 

Again, the value of £■ does not explicitely enter, except through the 
requirement that for the perturbation approach (known as the guiding center 
theory) to be valid the problem must satisfy "Alfven’s criterion" 

( S/ B) (l>» i /Ux ;j ) « 1 

An implicit 6 may be "made visible" by the following device. Consider 
a Hamiltonian with slow time dependence: one may artificially introduce £ 
into its time derivative by writing 

O^t ■ C OH/T>(6t) (6) 

Since 

^H/^et a 0 ( 1 ) 

this notation clearly displays the fact that the term is of order £ , 
and for this reason the Hamiltonian ( 2 ) is often written 


H 


H(p, q, £t) 


- 6 - 


A similar device may be used when there exists a slow dependence on 
dynamical variables; this can be quite useful in arranging the terms 
according to their orders in £ , but two things rust bo remembered. First, 
because of the way in which E is introduced, expressions of the k-th 

k 

order which have a factor £ standing in front of them, will also have 

-k 

"hidden inside" a factor £ . Secondly, because a definite value of £ 

is never stated, such factors must be cancelled out before the final 
result is obtained. 

An example may be useful here. Suppose a one-dimensional motion is given 
with a Hamiltonian that has a slow dependence on t , and it is also given 
that if this dependence is "frozen" (limit t - 0 ) the motion is periodic. 

The solution of such a motion usually begins with a canonical transformation 
to new variables (P, g) which are th<t action-angle variables of the unpertur- 
bed motion. If S is the generating function of this transformation, which 
in general is also slowly dependent on t , then the new Hamiltonian H' is 

H'(P, &) = H + 'Ds/'bt 

= H + £ "DsA (£ t) (7) 

H .(°) 4. £H’^ 

In the transformed Hamiltonian, the first order correction H'^^ has a 
factor t preceding it, but this factor is artificial and is balanced by a 
factor that is "hidden inside" the term, as is evident from the deri- 

vation. In practice, these factors must be cancelled before, say, the canonical 
equations of motion are used. 
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THE POINCARE - VON ZEIPEL METHOD FOR SLOW TIKE DEPENDENCE 

Consider a canonical system with 2K variables which has a slow dependence 
on time. We assume that the Hamiltonian H may be expanded in powers of £ 

H - !» ( Q ) 

We further assume that the Hamilton-Jacobi equation for £ 0 has been 

solved and that the transformation derived by it has already been applied, 
deriving as action-angle variables Cor the unperturbed motion 

( J » -A. ) = (p lf q x ) 

and giving 

= J t^/2 7\ (9) 

with (-0= u)(t) the slowly varying angular velocity. In the limit £ 0 , 

evidently, u) is a constant and so are all the canonical variables, except 
for which is then linear in time. 

To "solve" the motion we now seek a near-identity canonical transformation 
to new variables (P, £) , with 

( j\ n*) = (Pj. q 1 ) 

generated by 

<r(£,. q, t) = 21h li ♦ 2L £ k <r^(£» q. t) (w>) 

such that the new Hamiltonian H* does not depend on .TL • This is 

(15) 

somewhat similar to, but simpler than, an approach advocated by Gardner' 

( 2 ) 

and investigated by Contopoulos' , in which the same result is obtained by 


8 


a succession of canonical transformations, each of which pushes the elimina- 
tion of i”L* from H* one order higher. 

If H* is expanded in a manner similar to (8) and the time derivative 
is expressed as in (?)» one obtains 

22 t k h*M(p, a, t) » 22 £- k H ^(p> a> ♦ 

+ 22 e k ( TXT^^V-a fct ) (11) 

This equation contains 4N canonical variables, but half of them can be 
eliminated by means of the transformation equations 

Pi - p i + 22 £ -“ n>0' (B V'3q 1 ) ( 12 ) 

“i = q i * 22. £ " 0<r(m)/7) v (13) 

To facilitate the elimination it is best to follow a method introduced by 
Musen^ and use expansion operators ( * demotes operation, P 

etc. are gradient- type operators): 

H*^(P, a, t) = H*^ ( P, q + Jj.* t ) 

(14) 

= .*pjZ2 t"( '3<r ( “V'3EMV'3 3)} * H* (k) (p, q, t) 

■ EL=1 J 

= 22 £ m T< m) * K* (k) (P, t) 

m=0 “ 

where 


& - fmt 


a m 


i 


r 


i - 


9 - 


T<°> - 1 

T (l) - ('&<r (l) /‘3P 1 )(T)/'»q 1 ) 

l 

t* 2) * Hj ('SO' (2 Vt)P i )('^/'3q 1 ) + (15 

l 

♦ i 2Zl <^0' (l) /7>P i )(^<i'( 1 V'&P J )('7) 2 /T 1 qp qj ) 

i. i 


etc. Similarly 


H^(g, 

a. t) - £.“ 3 (m) . H (k) (p, q, t) 

fflnO 

(16) 

where 



s <°) = 

1 



£ (^<r (l) /W('aA>P i ) 

i 

(17) 

^ - 

7 1 , (^cr (2) /'jq i )0/7)P 1 ) + 



i 


* * £ ^ cU) /'3< i )('!)<r( 1 V7)<j J )(T) 2 /'aP 1 ^p J ) 

1J 

and ao forth. Substituting all thie in (ll) and collecting terms associated 

with £ gives 

k k 

*(■) . h* (I c -“) - £ ■ w ** (to * ) ♦ 

I0?O El*0 

( 18 ) 
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*(k) (k) 

The terms with m=0 simply equal H v ' and U v and will be taken 
outside the summation. The terms witi m*k also have simple form, for in 
general 

s (k) =» ( 0 )CT^ k V^q 1 )(^A> P i ) ♦ N^ k ^ (19) 

(k) 

where N contains only terms with at least two differentiations. 
Substituting ( 9 ) then gives 

s ( k ) # H (°) a ( w/ 27 \ ) (20) 


Because the transformation reduces to the identity transformation in the 
limit of vanishing £ , H*^ equals and due to ( 9 ) it satisfies 


T ( k ) # H *(o) ^ Q 


since operators involve only differentiation by the q^ , which 

does not contain. One then obtains the basic recursion relation 


(cO/2Jl ) Q<J-( k V<)A - H # ( k \p, q, t) 


with 


k-1 


(k) . r ( t ( b) . s (m) # H (k-m) 






A (k) (P. q, t) (21) 


^ - 'Jcr /-"a £ t 


( 22 ) 


depending only on orders lower than the k-th. If A. enters only as an angle 
variable with period unity, any function P(P, q, t) may be resolved into 
an "averaged" part 


* i ts 
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1 

? > -■ { F dfL (23) 

0 

and a "purely periodic" part with taro average 

(F)p. P - r - <»> 

Tha darivativa of a purely pariodic function ia also purely periodic 
and therefore, 30 -a 

'■fccr^V'Wl - Q/'Dil {(0^)^ ♦ (24) 

since tha contribution of vanishes. On tha other hand, 

does not depend on iT , so one gets the recursive relations 

H*^ - - <^A (k ^ (25) 

^Cr( k )/Qa » (2T\/A) ( /\^) (26) 

JWT 

Once these are solved, the calculation may be advanced to the next order, 
EXAMPLE 1 THE HARMOMI C OSCILLATOR ^ 

The Hamiltonian H" of a harmonic oscillator with a slowly time-dependent 
angular velocity c5(t) is 

H" = (l/2m) [p 2 ♦ ol 2 m 2 Q 2 J (27) 

If one "freezes" the time dependence, one can solve the Hamilton-Jacobi 
equation and derive a canonical transformation to action-angle variables 
(j, A ), generated by 


r““TT 


r 
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j ^(Ju>b/tt ) - B 2 Q 2 ^ 


dQ 


(28) 


Following this transformation, the new Hamiltonian H becomes 


H 


Jti/27T * £ J (oi'/4Tioi ) sinU^I-O.) 


(29) 


where the dash henceforth signifies the operation '7>A<) ( £ t) . Let (T of 
(lO) generate a transformation to (j*, IT*) such that all orders H # ^ 
of the new Hamiltonian are independent of IT * . This, combined with the 


fact that in the present case the only differentiation performed by T 


(m) 


of (19) is , allows all such operators to be ignored except 

«(o) 


for r 


A further simplification is obtained by noting that H contains only 
two orders, both linear in J 1 using the argument of (19) this gives, for 
che terms of (18) depending on H 


2, S (m) . 


3^), *. S^ k_1 ^ 


(30) 


• (^/2j\)'bcr^ k V'3n + ('a<r^ k " 1 ^xT.)(« , /43i(j) sin (4)in) 


In what follows, we will for conciseness write J instead of J , restoring 
the superscript — if necessary — only at the end. In analogy with ( 2l) 
we then obtain as the basic recursion relation, for k > 1 


( o)/2 1\ )(D CT V'Tn ) - H 


*(k) 


(31) 


- ('T>(r^ k '’ 1 V'T>rL)(c^'A7U^) sin(4^Sl) - (CT^- 1 ))' 


t *. 


TT 
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Using (18) directly for ks=0 , one simply gets the equality of and 

, while for k*l this yields 


< oi/2-xi ) i>c r^Vaa ~ h # ^ 


- J (ui*/47u3) sin ( 4 Jtn) (32) 


from which 


H* (l > « 0 


(33) 


<rM . 


j (cO’/earco 2 ) cos (4Jia) 


( 34 ) 


Higher orders, derived by the use of ( 31 ) t are 


H 1 


►( 2 ) 


(J/16JT) (^-)% 3 


cr 


( 2 ) 


p 2 l 

(J/647T) (<d'/ur) sin (Snn.) - >(35) 


- (j/l6T\u> ) (cJ'/u> 2 )’ sin (4flft) 


►(3) 


0 


CT ^ = - (j/384JT)(iW^ 2 ) 3 cos (I2hib) 

- (J/1283T w ) j^( Oi'/u) 2 ) 2 J cos (8 C\.) + (j/l28JT)(co'/‘A> 2 ) i cos (4Jlfl) 


- (j/327\^) J^( cO */t0 2 )*/oJ cos (4^n) 


( 36 ) 


Note that any term in an 0( £. ) expression cor tains the dash operator exactly 
k times, corresponding to the factor £ "hidden inside." 
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At this stage eqs. (12) and (13) could be used to express (J,.TL*) in 

terms of (J*, A ), up to O(C^). In fact, expressing the result in this 

manner, in terms of "mixed variables", is not too useful, and it pays to 
"invert" the result and express (j*, i\ # ) in terms of (J,iTL) , or vice 

versa. The shortest way to achieve this is by means of the direct transformation 

(17) 

technique' ' • If 

y a (p» q) 

are the "old" variables and 

z, » (P, 


are the "new" ones, and if the relation between the two sets has the"direct" 
form 


k = y + XU £ k $ (k ^(y) (37) 

k=l - ~ 


(k) 

then for this to be a canonical transformation, £ v must have the form 


(k) 


(38) 

where V is a gradient operator in "conjugate phase space" 


y s (q, -p) 

/(k) (k) 

the Tv are arbitrary functions and fr are prescribed expressions 

involving lower orders. In particular, if (37) is the "direct" form of 
the transformation generated by (10), one may choose 


^(,^(y) = - CT (k) (p, q) ■ - CT (k) (y) 


(39) 


,(k) 


(i.e. P is everywhere replaced by p ). The corresponding f'* 7 is 

k-1 


,(k) 


(“L n^r-(k-m) 


- T u'“'» 7 <r 

m=l 


(y) 


( 40 ) 


1 


T 


¥ r 


15 


with 

expansion 

operators depending only on the momentum-like com- 

ponents 7\ 

(m) of 





5 (o) 

a 

( s x , ... 

$m b) . 0. ... 0 ) 

with 

„d) 

3 

,(l) vr 
1\ * V 

(41) 


0 (2) 

a 

j) ^ 2 V ♦ 



and so on. 

Of particular Interest is the derivation of the adiabatic invariant 


z 1 = J* (42) 

which will now be outlined. 

To obtain G^\y) one simply uses the expressions (34 )— ( 36 ) 
without restoring the asterisk superscript (as was originally planned). 

To derive (4l), note that only one component of canonical momentum enters 


the calculation, so that 


o< l > - 

('"b/'D j) 


(43) 

u< 2 > = 

J # ^ (o)/dj) + i (j*^ 1 ^) 2 ^ 2 /^J 2 


The second-derivative terns may be safely ignored, since all orders of 
(k) 

O' ' used here are found to be linear in J . Finally, the components 
of the conjugate gradient V contributing to z x are simply 

/ ’iar (k) Ay 1 = 


^<r (k V'Dn. 
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For the first order, 



vanishes and one obtains 



0<r (1 V^a 


■ J (u)'/2 u> 2 ) sin (4^\a) 


( 44 ) 


The next terns are 


= (J/8) (d'/rt 2 ) 2 ♦ (j/4b>) (cO'/ui 2 )' cos (4J)a) 

(45) 

= (j/16) (cO’/tO 2 ) 5 sin (4T\a) - (J/8A) J( i^uS 2 )'/u}] sin(4^a) 


THE ,, QLD lt NOTION OF ADIABATIC INVARIANCE 

(«A 

In some texts of mechanics s r and in the older literature, the definition 
of adiabatic invariance differs somewhat from the one given here. The alter- 
native definition is usually applied to one-dimensional systems (though 
generalizations for several dimensions exist) and is as follows: 

"Given a slowly perturbed periodic motion, consider the action 
integral 

J =6 ® p dq (46) 

evaluated over one period of the unperturbed system. As the 
system is perturbed, an "instantaneous" J may be evaluated 
at any time by "freezing" slowly varying quantities. Then J 
has the property of adiabatic invariance: if the system undergoes 
a finite perturbation — e.g. , a finite change of the Hamiltonian 
from to Hg — the corresponding change in J may be made 
arbitrarily small by stretching out the perturbation over a 
sufficiently long time." 


The action variable J of (46) is the same as the zero -order action 
variable with which the previously-developed perturbation scheme begins , 
tut its "adiabatic invariance" differs in two respects from what was 
earlier defined as adiabatic invariance. First, there exists here no 
hierarchy of invariants each of which is conserved to some specified order 
and secondly, the definition concerns itself with the cumulative change 
in J over a long period in time. In fact, this property does not follow 
automatically from the definition of adiabatic invariance used amrlier. 

It is nevertheless an extremely useful property, since it allows one to 
derive, using only the unperturbed variables, a quantity with long-term 
invariance properties, without even specifying the perturbation . 

Since J is the zero-order part of J* , we may use (12) to obtain 
(compare also eq. 44) 



e'cXr^Vm + o( t 2 ) 


(47) 

» j + e j* ( i) + o(e 2 ) 

As in (12), <J^ means <T^(«J*, A , t) i since J* is a constant 

of the motion, only A and the slow direct dependence on t contribute 

*(l) 

to the variation of the first order correction J . The basic reason 
for the "long-term adiabatic invariance" of J, stated earlier, is that 
by the arguments of eq. (24) J*^ is purely periodic in A , and there- 

fore "nearly" purely periodic in t . Over long time intervals, its vari- 
ation is therefore bounded, causing the long-term conservation of J to 
be better than might otherwise be expected. 


16 - 


To demonstrate this, expand (47) to 

J a J* - t) ♦ 0 (e, 2 ) 

(48) 

= j* - £j*^(j\n., o) - e 2 t e.t) + .. + o(£. 2 ) 

Let a time T =* 0( £~*) pass. The first term on the light is conserved, 

while the second one will vary only through the variation of JPt. • Since 

the dependence of this term on Dl is periodic, the resulting contribution 

is bounded and due to the factor preceding it, of order £ . The next term is 

also 0( £ ) and the same holds for higher terms in tne expansion of the 

slow direct time dependence of . The 0( £ ) terms may contribute to 

dj/dt a term of form t S' » but its contribution to the total change 

of J will again be of order £ : 

T 

£ J" 4> dt ■ 0 (c 2 T) = 0( £■ ) 

0 

Hence the long-term variation of J is 0(£ ). 

The variation of other dynamical quantities, on the other hand, will 

be finite. For instance, for H 

T T 

AH = j (dH/dt) dt « J ('DH/'Dt) dt 

0 0 (49) 

» T^H/^aver. = (kT) (W^(et)) ayer 

and each factor here is 0(l). Thus by making £. arbitrarily small, but keeping 
T = o( C~^)( the variation of J may be made as small as is desired while 
that of H remains finite* 
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THE POINCARE - VON ZEIPEL METHOD FOR SLOW DEPENDENCE 

ON CANONICAL VARIABLES 

Let a perturbed periodic notion be given, represented by a Hamiltonian 

H - 22, t k 0 (50) 

with (p { , q,) the action-angle variables (j, ft) of the unperturbed motion) 
Since we have already derived methods dealing with slow time dependence, we 
will simplify matters by not including such a dependence here* The motion 
represented by alone is assumed to be periodic and soluble: we shall 

not require at this stage that has the form (9), but we note that it 

must be independent of -TL , since J is a constant of the unperturbed motion* 
Instead, we shall assume that the canonical variables y^ fall into 
two groups: "normal" variables for which 'li/T) y^ maintains the same order 
in S' and "slow" ones for which it raises the order by one level* It is 
useful to define parameters that distinguish between the two groups: let 
equal 0 or 1 depending on whether q^ is normal cr slow, and let £ ^ 
play the same role for . One can then define 


Qi " 

S' 9-^ 

(51) 

ip i - 

L Pi 

(52) 


so that (for example) OH/DQ i and OH TP ^ are always of the same 
order as H itself. 

As before, let a generating function 

<r (p, q) = YL p i 'k * L (53) 

1 1 

•Jt K* — 


1 


K 


J 


define a near-identity transformation to a new 


canonical set (P, , with the new Hamiltonian H* independent of the 

transformed angle variable i\* (,a tens with k m 0 could be included, 

but since it may not depend on ft it is not useful here). Again, the basic 
equation is 

H*(P. 3 .) - H(p, q) (54) 


and this again is expressed in powers of £ and expressed solely in terms 
of (P, q). Since p no longer appears, it helps to redefine 1P t as 

(P i =■ ^ (55) 

and this will be the definition used in the remainder of this section. In 
analogy with (14) one finds 

H* (k) (P, a) ** «P j H £“ Hi ('J>0' ( “ ) /'»P 1 )(^>A qj | . H* (k) (p, q) 

v 0=1 i J 


exp 


L 

i,m 


m+ S; + Vi 


(W“V'aP i )('a/aa i ) j * H #(k) (p, q) 


«P (■»<>- (m ' 5l ‘' rt) /'an 5 1 )('&/'3fi i ) . H*( k) (P, q) 

m i 


= JL £.* ,,(b) <* »* (1c) 

m=0 


(56) 


with V 


(m) 


suitable operators and 



vanishing for all non-positive 


values of m . Expanding the exponential gives 
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T (0) . l 

(57) 

- 22, 5i_ri V'Dfl* 1 )('J/' , *® 1 ) 

i 

and ao forth; because thasa oparators are axpraeaad aolely in terma of (P ^ 
and <Q ^ , thair action on any function maintains tha ordering in powara 
of £ • 

Similarly 

H (k) (p, ,) - £ t B « W * H (k) (£, q) (58) 

m=0 

with 

R^ » 1 

(59) 

r ( 1) = 22 ('‘>< r(1 ~^~^ Va<Q i )( / ^/^n , i ) 

i 

and ao forth. Substituting these operators and collacting terms associatad 
with £ ^ then gives, in analogy with (l8) 

22, { V (o) . H* (kH,) - R (b) . H (k -“ ) } - 0 (60) 

m*0 

Again, the terms with m=0 and mak are separated. For the latter terms 
one gets, in analogy with ( 19 ) 




22 - 


z> 


(k-5*- Yi) 


1 V"<>(P i )('c>/ > c><Q i ) ■*■ M' 


(61) 


* L ^^ (k “ S -^ ) /^Q i )(W7)IP i ) ♦ 


with and N^ involving only lower orders. By taking k=0 in (60) 


one again fiiids 


E *(°) , H (0) 


(62) 


so that (60) becomes, for the generfl case 


H* (k ) + 


ZZi [(-7><r (k-?1 “ 5i V'dF i )(W°V'3© 1 ) - 

- ca<r (k " Si " <i) A>© i )('^H (0) /T)n > i ) | 


(63) 


where 


0< k > . H^ k ^ ♦ 


^ [r U) . H (kH,,) - V W . H* (kH,) } ♦ (M (k) - N (k) ). H ( 


(64) 


involves only given functions and lower orders. For every k a relation 
of this type is obtained, constituting a k-th order recursion formula for 
the derivation of and <J~ . 

Now the action-angle variables associated with the zero-order periodicity 
(and used here in "mixed" form) 

(A n. ) a (F r q x ) 


i . ft 


23 


are assumed to be "normal", so that the left-hand side of (63) will include 
a term 

j*) 

(k) 0f 

If 0"' appears nowhere else, the equation assumes the formal) and 

is solved in the same manner. 

On the other hand, if appears anywhere else in (63), it may 

not be possible to derive it, for the equation then becomes a partial 

(k) 

differential equation for O’ ' ' . To prevent this from happening, it 
is required that for all (p^, q^) appearing in other than the 

action-angle pair 

Si + Yi > 1 (65) 


Hence the recursion can be carried out if: 

(1) has normal dependence on p^ but does not depend on q^ 

(2) H^ may depend on any "slow" variable. 

(3) H^ may depend on any "normal" variable, provided its 
canonical conjugate is "slow". 

Furthermore, it may be shown by extending the present calculation that 

/ \ (-l) -1 

(4) H may include a term H' ' of order g. , provided it depends 
only on slow variables having slow conjugates. Such terms are the*' 
transformed intact to the new Hamiltonian. 


As an example, the Hamiltonian of a charged particle in a time-independent 


electromagnetic field , in the regime of guiding-center motion, may 

(i5)(a$) 

be brought to the form 
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H 


p 2 2 /2m ♦ p^tO /2JT + e _1 a<^ 0 ^ + P(£) 


( 66 ) 


Here (p^, q^) are canonical variables associated with the rapid gyration* 

(p 2 , q 2 ) represent the motion along field lines and (p^, q^) describe 

the identity of the guiding field line* which changes slowly with time 

(in the references* subscripts 1 and 3 have reversed meanings); the 

variables (p 1# q^, p 2 ) are "normal" whereas the remaining ones are "slow". 

Furthemore* the gyration frequency co and the lowest order <f> ^ of 

the electric potential are both functions of slow variables only. The last 

term is of order £ -1 , since its derivatives are proportional to the 

of the 

components of the lowest order^electric field B, which are of order 


unity. 

Evidently H meets all the previously stated conditions except for ones 
if <f > ^ contains q 2 * condition ( 4 ) is violated, since p 2 is not 

slow. One therefore must impose an additional requirement that q 2 

vanishes: this reduces to the well-known restriction in guiding center 
theory that the electric field may have no zero-order component parallel 
to the magnetic field. 


DIRECT canonical transformations with slow variables 

(k) 

The generating function CT gives the transformation equations 
as in (l2)-(l3), in "mixed" fox®. To bring them to the "direct" form (37) 
it is useful to generalize (38) for cases in which slow variables are 
present* and this will now be done. 




1 fi. 


1 — rt 
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Let 


y * 



q) 


be a canonical set and 

y = (q. -p) 

be its conjugate^). One may now define an index equaling 0 or 1 

depending on whether y^ is "normal" or "slow” , and an index A ^ which 
has a similar relation to y^ . With this notation it is poesible to define 
vectors Y and Y satisfying relations similar to (5l) and (52) 


= £ r ‘ y A (67) 

\ y ± ( 68 ) 


As with quantities defined in (5l) and (52), Y^ and ^/DY^ are 

always 0(l). 

We now seek the condition for a near-identity transformation (37) to 
be canonical. Actually, in what follows the recursion may still be carried out 
even if the transformation is not a near-identity one and £ terms are 
included (satisfying appropriate conditions) but we shall not develop this 
possibility here. One then finds, as a condition for canonical behavior 
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Expressing derivatives in tezms of Y and Y gives 


[a, b] =* y Q ) 

s 

= Yh (Ha/-T>Y s )(^!b/^Y a ) 

8 

In particular 

[a, = e Al (T>aA\) 

Thus 


( 70 ) 


( 71 ) 


0 ■ Ef^’ ( ^ k) /*V - £* 4i 

lc*l 

( 72 ) 

♦ 22 '£ k+As+r ‘ 2D (^“Vav^j^^v} 

s m=l 


Dividing by £^*+ ^ 
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' i $[ k) /7>?,) - e k ‘ 4i (^S ■ k V'3? 1 ) + 

A. A A _A_1 I X j 1 


*■ A i+ A^L 


(73) 


k-1 a 

2^ £ k + r, + A s -4-Aj J-’ > 

8 IQsl 


It is useful at this point to redefine k for each term so that all powers 
of £ become g end also to replace m by 

M » m - A L (74) 

Because the exponent of £ differs for each term in (73) » the new summation 
over k will begin at a different value for each term; this summation limit 
may, however, be uniformly set equal to 1 if it is assumed that £ ^ 
vanishes for non-positive values of u . With these changes, (73) gives 

o = Y s £ \ Vmj) - + 

k=l l 

(75) 

k- A i ~l 

♦ 2Z JZ (^< M+A ' J /T>Y 8 )(T>^ k - M - rt - is+ ^ ) /T>Y 8 ) • 

M»l- A ^ a j 


This suggests the introduction of new "staggered" vectors 



$ 


(k+&l) 


(76) 


r 


* t 
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■2 oo 




(k) 


if - 0 


(77) 




(k) 


5 £ 


(k+l) 


if - 1 


It is also ussful (in analogy to what was dons ' n ref. 17 ) to introduos 
a curl operator in 7 space. With this notation (75) nay be rewritten as 


< * T * t\i 


(78) 



k-l-r - 4 *4. 

8 8 J 



8 


M=.l 


('2)^ M) /'ai 8 )(^^ k ‘* - rs_A$) /'ai s ) 


,(k) 


This equation may in principle be used to derive recursively, 

but this turns out to be a rather inconvenient approach. It is more 

useful in determining the degree of arbitrariness associated with a 

near-identity canonical transformation of the font (37). Let two such 

(m) 

transformations be given, characterized by "staggered” vectors ' 

and x (m) w’.aich are identical for orders up to and including the (k-l). 
For the k-th order one finds that the right-hand side of (78), which 
depends only on lower orders, is identical for both expansions, giving 


V y X 


( ^ - * (k) ) 


(79) 


TT 
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from which 


. ,W * V Y X (k) (80) 

Thus the arbitrariness in specifying the canonical transformation at each 

(k) 

level of ^ ' is contained in the gradient in Y space of an arbitrary 

(k) (k) 

scalar . The general form of y ' for a canonical transformation 

may be written* in analogy with (38) 




(81) 


vh.r. p(^) is a vector Involving orders of lover than the k-th 
and constitutes one particular solution of (78). In the following sections 
two such particular solutions will be derived, analogous to those found in 
ref. (17) for small perturbations. 


DERIVATION BASED ON (T(P, g) 


Let a near-identity transformation of n » 2N variables 

y = (p, q) z = (P, £) (32) 

be given by (37), and let a generating function (33) be assumed to produce 
the same transformations via the equations (12) and (13). In what follows 
the relation between (37) and (33) will be established in a way resembling 
what was done in ref. (17) for the case when no slow variables are present. 
As before, the calculation may be broadened somewhat beyond what is done 
here, since the method only requires that canonical momenta transform in 


near-identical fashion. 
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With the notation of (5l) and (55) eqs. (12) and (lj) giva 

p 1 - p t - £ k (^> cr^V-a q^) 

• p i (/ *° W/ v 

\ - ^ v 

■ qi + 2Z ektri(1><r(k)/ ' 1>1Pi) 


( 83 ) 


(84) 


All functions on the right depend on nixed variables (P, (j) ; to introduce 
a dependence on y , it is useful to define "partial vectors" adding up to 


( ... • °> — 0) 
( 0, ... 0, J^ k) . ... $* k) ) 


(85) 


Prom this 


p i + /\ f kjT i k ^ 

■k *21 & k0 *A ( ? ) 


(86) 


(87) 


If the vectors Jr are known* they may also be used to expand any 
function of nixed variables (P, q ) in terms of y * e.g. 


1 
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r (£, q ) - T(j l k J \ (k ') 

* trp (y -1 . t k ( JT ^ 'b/* y)^ • J(jO 

■ “P fc £. k 22 ( 5 (k_ r ‘ T>/t» V j « *(y) 

- 2Zi £- k t(k) * F( ? ) (ae) 

kaO 


vhin, if we collectively denote all "alow" components of y by £ and 

(k) 


all "normal" onea by £ and if 


31 i k) and J\ 

— m •“ r 

(k) 


composed of the corresponding components of T\ 
L<°> - 1 


denote veetorm 


a) 


J} ,* '<>/■»£ 


(89) 


Tip 2 } ‘i/'t r + D/DI& 

* * - r^ 1 C^/'^xH'JAx) 


aince ^ (i) (?) 

ami so forth. Note that £ is implicit, jv ' and T1 ' ' should 

-1 -2 

have factors £ and £ "hidden inside", since they are teamed up 

( 2 ) 

with the corresponding positive powers in (66). For the same reason L v ' 

-2 

should contain a factor £ and indeed, inspection of the last equality 

fit 

in (89) shows that all terms there have such factor . 

\ 


Substitution in (83) yields 
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\ • p A - y -1 t“ i* <o) • ('a<r (k v^<Q i ) 

' 1 " m 


P< * 


jT n ( e k f (T><r (k ' ,ri V - 9Q 1 ) ♦ 

V 

1c- ^ ^ 

+ £ 1 l(b) * ('J<r (k - B - r 

m=xl w/ 


(90) 


where all terae of CT are viewed as functions of y , i.e. with p 
replacing P wherever the latter originally appeared. This should be 
identical to (83) and therefore 


(k) 

JTi 


- - 2Zi itm> * 


/T)<Q 1 


(91) 

The highest order of L 10 appearing on the right is k-1 ( (J ^ only 

appears if $ 1 , for if It depends on "normal" variables, the trans- 

formation is no longer one of near-identity: this is the reason for the 
change in summation limit) and this is therefore also the highest order 
of JT ^ appearing on the right. Thus (9l) is a usable recursion relation 
for deriving J1 ^ . 

Expanding (64) in a similar manner gives 

k-1 

e£ k) - * (^CT^ k “ m " 


where the definition of IP ± reverts to (52) 


i I_ 


I 
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Now if 


then 


and if 


then 



Inspection then shove that (9l) and (92) may be combined to one equation 


i - y ± - Y e k | Dcr (k “ ^ 1 

^ — ‘ L m=*i 

(93) 


(k) 

Thn» dependence on & ^ may be removed by introducing defined 

in (76). Then, using the gradient operator in Y space, (93) becomes 

(k) » - V v CT (k) - ri W .?,(T W (94) 


v 


L 

mail 


Since it has already been established in (81 ) that is arbitrary 

within some gradient in Y space, the summation term represents a parti- 
cular solution of (78), 
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LIB TRANSFORMS Wlffl SLOW VARIABLES 


(16) 


If Ly is the operator denoting Poisson bracketing with a function W of 


the canonical variables 


Iv ( f ) “ 


(95) 


Then it may be shown 
y « (p, q) to 


( 17 )( 19 )( 20 )( 21 ) , 


hat the transformation from 


exp ( £ Ly) * y 


(96) 


(with the exponential operator defined by its series expansion) is canonical. 
In what follows the form of the Lie transform in the presence of slow vari- 
ables will be derived, again following closely the derivation for the 

( 17 ) 

simpler case when all y^ vary on the same scale . Let W be expandable 


in £ 




(97) 


and l9t op”®*™ 1 v ** ietu,ei throu « h 




( 96 ) 


If one. 


^defines (compare ref. 17, eq. 35 ) 


- W 


(k-l) 


(99) 


1 IB _ I 
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then 


>* = - L H '* (k) /'» y a )( '<>/'» y B ) 


k ZL ("bX (k " r ’" As ‘/^ s )(' 4 AY 8 ) 


where the lower limit of k in the last summation may be chosen as 1 
if quantities with negative or zero index are understood to be zero. This 


gives 


L (iox) 


Expanding a typical component of (9C) 


■, - {> * ‘E 1 " ■v"* * *£* k V k> > ! --} • % 
= N(k;) * ( c -” 1 " 1 v 


y & k M (k+ n k * i 


(102) 


where the M v all have the fonn 


L^(k) + H (k) 


(103) 


with N' K ' some operator involving lower orders. One then gets 


a * - _ « 
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r 

Ezr a 


Yt 


00 


$ 


(k+ &i) 


=t M^ k * r i. ^ * y. 


(k+r l + A L ) + M ( k+|- t + Ai) # y 


T>X k) /i)\ + n 


(fc+ r t + At) 


* Y, 


(104) 


which again is the sum of a gradient in Y space and an expression 
involving lower orders which (presumably) is a particular solution of (78). 

THE KRYLOV - BOGOLIUBOV - KRUSKAL METHOD WITH SLOW VARIABLES 


Krylov and Bogoliubov^^”^^^ investigated the solution of a set of 


n equat. 4 / ' vectorially represented by 

dy/dt » £ k g^(y) 

k=0 


(105) 


kith 


.( 0 ) 


(0, 0 , ... 0. g^ 0) ) 


(106) 


ensuring: that in the "unperturbed" limit £. 0 , y alone varies and 

all other components of y (to be collectively denoted by y ) are 
constant* It is further assumed that the unperturbed system is periodic 
and that y n is an angle variable appearing only in the angle-argument of 
periodic functions. The zero-order growth of y n is then assumed to be 
linear, from which follows that g^ may depend on jr but not on y Q • 


It 
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To eliminate the periodicity from this motion, Krylov and Bogoliubov 
used a near-identity transformation to new variables z , given in a direct 
form inverse to that of (37) 

y = ' z ♦ y. e k ^ (k) (i) (107) 

The new variables, which can be derived by a suitable recursive method, 

have the property that the equations by which they evolve do not contain 

the transformed angle variable z on the right-hand side but have the 

n 

form 

dz/dt =* £ k ^ k) (D (l<») 

The first (n-l) equations of this set, representing dj^/dt, then 
form an autonomous set not involving z^ and can be solved independently. 

If y represents a perturbed periodic canonical system with a 
Hamiltonian of the form (50), then the canonical equations of motion 
have the form (105) and the Krylov-Bogoliubov method o«n be used to 
eliminate the angle variable = fl . Unfortunately, unless precautions 
are taken^ 22 ^ 2 ^ , the z variables will in general not be canonical, so 
that the transformed variable corresponding to the canonical conjugate of 
y n will in general not be a constant of the perturbed motion, as is 
automatically achieved by the Poincare-Von Zeipel method. 

On the other hand, the Krylov-Bogoliubov method has a much wider validity 
and can be used in non-hamiltonian systems. A similar elimination procedure 
which derives the transformation in the form (37) has been devised by 


1 
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Kruskal' 0 '' , who followed it by the derivation, (for canonical systems 


only) of a constant J of the motion, obtained by an ingenious application 
of integral invariants (it is the same constant as is obtained by the 
Poincare-Von Zeipel method). 

Here the Krylov-Bogoliubov method will be generalized for the case when 
slow variables are present. As with the Poincare-Von Zeipel method, this 


allow the restrictions on the form of the zero-order equations — embodied 


( 0 ) »«• 
in g — to be eased. Specifically,, variables other than y are now 
— An 


allowed to have a zero-order variation and this variation (as in the 
canonical method) is passed intact to the ''reduced" equations involving z . 
The calculation will be done for the transformation (107) ; the treatment 
of Krruskal 's method, using (37), follows practically identical steps and 
will therefore be omitted. 


Following the notation of (89), let R and r denote the slow and normal 

(k) (k) 

components of js , and let <y tji and be corresponding components 

of 71^ . Substituting ( 107) in the left-hand side of (105) gives, with 
the definitions (67) (68) extended to z variables 


dy/dt 


dz/dt ♦ (H)^ k /'t>z B )(dz s /dt) 

k, s 


d*/dt + 2Z £* (W k ~ r * V'OZ^) J’ie." h 

k, s 


(■) 

8 


BIpO 


k-1 


dte/dt + ^Zi ^ Hj ( ^V k " " ~ ri V-ZJ 2 a ) h 


(») 


s 


( 109 ) 


maO 


Expressing a typical term of the right-hand side of (105) in terms of z, 
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gives , in a manner similar to (88) 


S 


(k), 


(y) - § (k) ( *+.Ce"l W ) 


m^U) 


° XP | Y2j ^ ^ i* 1 z s^} * 

1 m f s 

exp \J2jt-™ ( ^•'V' 5 Z B ) ] * « W (e) 


= £ t“ k W * g«(*) 


(no) 


The operators resemble those of (89) but with -v^ m ^ everywhere 

/ \ 

replacing Jr ' . Substituting preceding results in (105) gives 

r k-1 


^ — 1 


dz/dt = 

k=0 


2^{k w . « (k -° ) - rs)/r3z s ),h i o) } 


m=l 


s 


♦ (K< k >. g (0) ) 


£ * W g<°> 


a 


+ ^ ^kCp 5 


(k) 


(111) 


where the factor preceding the last term denotes that it be omitted for k = 0 
(in that case it is already counted as the terra involving K^) and where in 
the summation preceding this term h^ has been replaced by , which 

equals it since in the limit eqs. (105) and (108) coincide. 

Comparison with (lC8) shows that the expression in the square brackets 


i * 


« a. 


rt 


equals h^(jz) » and this equality forma the basis of the recursive deri- 
vation of and . 

The situation now resembles that of (63) : in order that the recursion be 

at all possible, unknown components of must appear in (ill) only 

once , otherwise the result is a partial differential equation and cannot 

(k) 

be easily integrated. One term which always contains ' is contributed 

by the last summation in (ill) and equals 

g^ 0) 

since gj^ does not vanish and z r , the transformed angle variable, is 
M noxmal". No other appearance is permitted, hence 


4 0) - 0 


for 


r =0 

8 


( 112 ) 


or, stated in words, "Only slow variables and the principal angle variable 

are allowed to have a zero-order variation." 

/. \ the term containing 

In addition, * could enter through^ K', which has the form 

(compare eq. 89) 




(l 


GO. 


O/'Ds 


N< k '' J « g (0) 


( 113 ) 


(k) 

with N v ' containing lower orders. No problem arises here provided 'vj ^ 

(k) *" 

is derived first and ^ ^ only afterwards: because of (112), this term 

is absent in the first part of the derivation, while in the second part those 

^(k) 


components of 


that appear in it are already known. In either case one 


f 
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g< 0) ♦ h w (z) - > 


(k), 


(k) 


:u4) 


where (k > 0) 

v(k) . 


£ { K (m) . ^ ■ r ‘ V**.) h<“> 


m=l 


8 


(115) 


♦ (K< k W 0) ) + g W 


ZL (^^ k " rs) /T)2 s ) «i 0) 

8 ^ n 


depends only on lower orders. The Solving of (114) then resembles that 
of (21) 


CONCLUSION 


In the preceding sections the main methods of classical oerturbation 
theory have been extended to slowly (or adiabatically) perturbed systems. 
At the same time, the basic concepts associated with such systems (e.g. 
adiabatic invariance and implicit 6 ) were examined and clarified. 

The restrictions on the forms of the zero-order equations for slow 
perturbations have been derived and are generally less severe than for 
small perturbations. The extension of the methods themselves is relatively 
straightforward, involving mainly the shifting of indices for quantities 
corresponding either to slow variables ot (as the case of to 

variables with slow conjugates. With the use of expansion operators the 
treatment is only slightly more complicated than for small perturbations. 


i 


K 


f 
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